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A  CANCELLATION  TECHNIQUE  FOR  REDUCING  BACKGROUND 
NOISE  WITHIN  TURBULENT  FLOW  ENVIRONMENTS 
CHARACTERIZED  BY  PIPES  AND  ANNULI 


INTRODUCTION 


Wall  pressure  fluctuations  generated  within  turbulent  fluid  flows  axe  of  practical  as 
well  as  fundamental  significance.  Besides  adding  to  our  interpretation  and  understanding 
of  the  physical  phenomena  present,  these  fluctuations  can  couple  with  various  structures 
containing  the  flow,  driving  them  in  such  a  way  as  to  cause  structual  failure.  Depending 
on  the  design  of  the  experimental  facility  and  the  properties  of  the  flowing  fluid,  the 
background  sound  level  can  range  from  small  to  very  large  amplitudes  when  compared  to 
the  turbulent  pressure  fluctuations  of  interest.  Hence,  discriminating  in  some  way  between 
the  two  levels,  better  enables  the  researcher  to  reduce  the  contaminating  background  noise 
and  accurately  measure  the  pressure  fluctuations  of  interest. 

There  have  been  a  number  of  approaches  in  the  past  in  investigating  ways  of  reducing 
the  effects  of  facility-generated  noise  from  interferring  with  turbulent  flow  experiments. 
The  simplest  method  is  to  design  the  facility  with  sound  isolators,  absorbers  and  linings 
such  that  extraneous  noise  is  minimized  as  much  as  possible.  A  number  of  investigators 
(1-7)  have  reported  on  experiments  utilizing  such  techniques.  This  has  resulted  in  a  fair 
degree  of  success  in  reducing  the  environmental  noise  within  tolerable  limits,  usually  on 
the  order  of  the  turbulent  fluctuations  being  measured. 

Another  technique  was  introduced  by  Willmarth  (8),  who  recognized  that  the  typical 
experimental  facility  contains  local  source  mechanisms  of  acoustic  sound  such  as  vanes,  dif¬ 
fusers,  bends  in  piping  and  pumps.  Willmarth  showed  that  a  finite  correlation  was  obtained 
for  negative  time  delay  between  flush-mounted  transducers  separated  in  the  streamwise  di¬ 
rection.  By  combining  this  temporal  quantity  with  the  spatial  separation  between  the  two 
sensors,  a  velocity  was  derived  equal  to  the  sound  speed  (relative  to  the  flow  velocity)  in  the 
fluid.  Willmarth  interpreted  this  quantity  as  flow-generated  acoustic  noise  and  subtracted 
it  from  the  original  output  to  eliminate  its  contribution  to  the  overall  signal.  However,  ap¬ 
plication  of  this  method  has  the  disadvantages  of  being  both  tedious  as  well  as  introducing 
an  additional  source  of  error  by  massaging  (correcting?)  all  data  prior  to  any  statistical 
analysis. 


MuuKript  improved  September  19,  1988. 


1 


A  third  technique  employs  a  signal  processing  scheme  introduced  by  Wambsganss 
et  ai.  (9).  and  later  refined  by  Wilson  et  al.  (10).  Utilizing  three  pressure  transduc¬ 
ers  located  in  a  plane  perpendicular  to  the  flow  direction,  this  method  employs  simple 
subtraction  in  order  to  separate  the  signals  due  to  the  turbulent  boundary  layer  pressure 
fluctuations,  background  acoustic  and  electronic  noise  respectively.  Both  Wambsganss  and 
Wilson  utilized  this  method  to  study  turbulent  flows  through  annuli. 

This  paper  reports  some  refinements  to  this  signal  processing  technique  which  mini¬ 
mizes  the  error  in  the  resulting  estimates  of  the  wall  pressure  fluctuations.  Experimental 
results  are  presented  for  pipe  and  annular  flow  geometries  showing  the  ease  with  which  the 
method  may  be  applied. 


CONCEPTS 


We  assume  first  that  the  turbulent  boundary  layer  wall  pressure  fluctuations  are  sta¬ 
tionary  ergodic  processes  (11).  We  further  assume  that  the  instantaneous  signal  from  each 
flush-mounted  pressure  transducer,  F(x,t),  can  be  represented  as  the  linear  superposition 
of  three  components: 


Fi(x,t)  =  Ti(x,t )  +  Ai(x,t)  +  E,(t).  (1) 

Here,  T,  A  and  E  represent  the  contribution  to  the  total  signal  due  to  the  turbulent 
boundary  layer  induced  pressure  fluctuations,  the  facility  generated  background  noise  and 
electronic  noise  respectively.  As  indicated,  the  acoustic  and  turbulence  signals  are  functions 
of  both  space  and  time,  whereas  the  electronic  signal  would  be  specific  to  the  transducer 
itself  and  result  in  time  dependence  only.  Using  standard  definitions  for  the  correlation 
between  any  two  time  dependent  signals  over  a  finite  time  P,  we  have  the  following  relation: 

1  fP 

Ri^FiFj  =  lim  —  /  F,{t)F}(t  +  r)dt.  (2) 

n  J0 

The  auto-correlation  is  defined  for  i  =  j  and  the  cross-correlation  is  for  i  ^  j.  If 
the  contaminating  noise  of  an  acoustic  nature  is  generated  several  hydraulic  diameters 
up  or  downstream  of  the  test  section,  then  the  resultant  pressure  can  be  assumed  to 
propagate  in  a  plane  wave  (or  axisymmetric  mode  for  pipes).  For  low  frequencies  and  stiff 
constraining  supports,  this  propagating  wave  may  be  considered  constant  in  magnitude  and 
phase  throughout  any  cross-sectional  plane  normal  to  the  flow  vector.  This  fact  provides  a 
means  to  conveniently  eliminate  this  component  from  the  instantaneous  transducer  output 
by  subtracting  signals  from  two  transducers  lying  in  this  co-axial  plane.  Referring  back  to 
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Eq.  (1),  we  see  that  the  other  major  component  is  the  turbulent  term  T,(x,t).  In  order  not 
to  eliminate  any  of  this  information  from  the  resultant  signal,  the  two  transducers  must  be 
spatially  separated  far  enough  to  assure  zero  correlation  between  the  boundary  layer  wall 
pressure  signatures.  It  is  generally  accepted  that  the  coherence  length  between  turbulent 
signals  in  a  transverse  direction  from  the  flow  would  be  somewhat  less  than  the  boundary 
layer  thickness  6,  hence  determining  a  minimum  length  scale  required.  When  the  above 
conditions  hold,  we  may  make  the  following  assumptions: 

(1)  Ai(t)  is  uncorrelated  with  the  turbulent  fluctuations  and  the  electronic  noise. 

i.e.,  (a)  /p  Ti(t)Aj(t  +  r)dt  =  0,  for  all  i,j 

(b)  /0P  Ei(t)Aj(t  +  r)dt  =  0,  for  all  i,j 

(2)  The  turbulent  fluctuations,  Ti(t),  are  uncorrelated  with  one  another  and  with  the 
electronic  noise. 


i.e.,  (a)  Jp  T,(t)Tj(t  +  r)dt  =  0,  for  all  i,  j 

(b)  /p  Ti(t)Ej(t  +  r)dt  =  0,  for  all  t,  j 

(3)  The  instrument  noise  in  each  transducer  is  uncorrelated  with  each  other, 
i.e.,  Jq  Ei(t)Ej(t  +  T)dt  =  0,  for  till  i  ^  j. 

When  these  assumptions  are  valid,  the  subtraction  of  two  signal  outputs  eliminates  the 
acoustic  noise  component  directly, 

F  =  Fa  —  Fb  —  Ta  —  Tb  +  Ea  —  Eg-  (3) 

At  this  point,  Wambsganss  (who  investigated  the  flow  in  an  annulus),  assumed  that 
the  electronic  noise  was  negligible  and  the  inner  and  outer  wall  pressure  signals  were  equal. 
Auto-correlating  the  signal  from  Eq.  (3)  results  in  the  following, 

1  rp 

R{t)ff  =  lim  —  /  [TA(t) -TB(t)][TA{t  +  t)  -  Tg{t  +  r)]dt  (4) 

P — 00  Jr  J  o 

=  TL&)  +  Tn&)  (5) 


=  2^(r),  (6) 
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where  the  overbax  refers  to  time/ensemble  averages.  From  our  understanding  of  turbulent 
flow  within  annuli  and  as  represented  in  experimental  results  (10,12),  Eq.  (6)  has  been 
shown  not  to  be  valid.  The  outer  wall  rms  pressure  fluctuations  exceed  those  on  the  inner 
wall  with  the  difference  increasing  as  the  ratio  of  outer  to  inner  radii  increases. 


Wilson  refined  this  method  by  first  auto-correlating  Eq.  (3)  while  maintaining  the 
integrity  of  each  respective  component. 


R(t)ff  =  TIa{t)+TIb(t)  +  E\a(t)  +  EIb(t) 


By  combining  Eq.  (7)  with  two  similar  equations  obtained  through  subtraction  of  the 
other  two  transducer  pairs,  a  set  of  three  equations  with  three  unknowns  may  be  obtained. 
By  solving  these  equations  simultaneously,  the  three  turbulent  terms.  T\A.  T'Bg  .  and 
T^q  could  be  obtained  once  the  zero  flow  signals  E\A,  E2BB  anc^  E&c  were  recorded. 
There  is  one  major  objection  to  the  use  of  this  method.  There  is  a  distinct  disadvantage 
when  the  turbulent  signal  is  of  the  same  order  of  magnitude  as  the  acoustic  signal  A(t). 
The  resulting  solution  requires  substractions  between  relatively  large  numbers  to  get  a 
smaller  one,  contributing  more  uncertainty  in  the  results.  However,  it  represents  a  large 
improvement  over  the  method  of  Wambsganss. 


The  authors  propose  an  improvement  to  this  subtraction  technique  which  involves  sim¬ 
ply  cross-correlating  two  different  subtracted  signals  from  three  coplanar  pressure  sensors 
as  illustrated  below. 


Fi  =  Fa  —  Fb  =  Ta  —  Tb  +  Ea  —  Eb  <  S) 

Fi  =  Fa  —  Fc  =  Ta  —  Tc  4-  Ea  —  Ec  !  9 ! 

1  fP 

R(t)f1f}  =  Jim  —  /  Fi(t)Fo(t  +  r)clt.  (10) 

P-oc  PJ0 

Assuming  the  assumptions  outlined  previously  remain  valid,  the  following  result  is  ob¬ 
tained, 


=  R(t)aa  =  T2a(t)  +  E2aa{t).  (11) 

All  that  remains  is  to  measure  the  auto-correlation  of  transducer  "A  ”  under  zero  flow  con¬ 
ditions  to  obtain  E2aa{t )  and  subtract  it  from  Eq.  (11)  to  obtain  the  turbulent  pressure 
auto-correlation  at  location  A.  One  may  note,  that  if  the  electronic  signal  EA(t)  is  suffi¬ 
ciently  lower  (say  30-40  dB)  than  the  fluctuating  turbulent  pressure  signal,  then  Eq.  ( 11 ) 
gives  the  desired  results  directly.  Of  course,  similar  subtracted  signals  may  be  derived  for 
the  other  pressure  correlations  R{t)bb  and  R(t)cc- 
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This  method  is  immediately  adaptable  to  either  longitudinal  or  circumferential  cross- 
correlations  to  ascertain  eddy  length  scales  and  convection  properties.  Referring  to  Fig. 
1,  the  following  relations  may  be  obtained  for  both  indicated  geometries  as, 

Fi  =  Fb-Fa  (12) 

1  rp 

R(t)FiFd  =  Jim  —  Fi(t)Fo(t  +  r)dt  (13) 

P~  °°  P  J0 

r(t)bd  =  nj7).  (U) 

As  indicated,  the  cross-correlation  of  two  partially  coherent  turbulent  signals  in  the 
flow  may  be  obtained  directly  with  no  additional  subtractions  required.  In  order  to  obtain 
the  spectral  content  of  the  turbulent  signals,  one  only  needs  to  Fourier  transform  the 
respective  correlations. 


EXPERIMENTAL  APPARATUS  AND  RESULTS 


The  signal  processing  technique  outlined  above  was  employed  to  measure  the  turbulent 
boundary  layer  and  facility  generated  noise  induced  wall  pressure  fluctuations  in  fully  de¬ 
veloped  pipe  and  annular  (external  wall  only)  flow  geometries.  Figure  1  shows  the  coplanar 
arrangement  of  three  transducers  located  in  the  external  walls  for  the  piping  arrangements. 
For  a  rectangular  channel,  one  would  place  the  three  transducers  as  represented  by  A.  B 
and  C  in  Fig.  1  in  the  external  wall,  side  by  side  in  a  transverse  direction  to  the  flow.  The 
flow  facility  is  a  blow-down  water  tunnel  operated  by  the  Laboratory  of  Computational 
Physics  and  Fluid  Dynamics  at  the  Naval  Research  Laboratory  and  is  illustrated  in  Fig. 
2.  There  are  two  test  sections  for  the  facility.  The  first  is  a  piping  arrangement  (used  in 
this  study)  and  the  second  (not  shown)  is  rectangular  in  cross-section  with  an  aspect  ratio 
of  18:1.  The  damped  pipe  test  section,  indicated  in  Fig.  2,  consists  of  1.S3  m  long,  19.1 
cm  inside  diameter  acrylic  or  polyvinyl-chloride  (PVC)  pipe  sections  supported  every  1.S3 
m  on  solid  cinder-block/concrete  foundations.  For  the  annular  configuration,  the  inner 
cylinder  was  constructed  of  a  3.81  cm  OD  brass  bar  which  was  supported  every  1.S3  m 
by  three  0.32  cm  OD  rods.  To  insure  that  a  fully  developed  turbulent  flow  was  present  in 
both  the  pipe  and  annular  test  sections,  the  pressure  sensors  were  placed  within  the  outside 
walls  more  than  100  hydraulic  diameters  downstream  from  the  ball  valve  (indicated  in  Fig. 
2).  More  detailed  information  on  the  facility  and  its  operating  characteristics  are  provided 
elsewhere  (13,14,15). 
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The  flush-mounted  sensors  are  capacitive/pre-amp  combined  units  (0.23  cm  in  diam¬ 
eter)  manufactured  by  Clevite  Corporation  of  Bedford,  Ohio  and  sire  similar  to  those  used 
by  Bakewell  et  al.  (2).  Each  combined  unit  was  matched  in  phase  and  exhibited  a  flat 
frequency  response  from  DC  to  10,000  Hz.  Each  experimental  run  was  limited  in  duration 
from  1  to  3  minutes  by  the  volume  of  the  upstream  liquid  reservoir.  Each  sensor  output 
was  recorded  on  a  14  channel  Hewlett  Packard  magnetic  tape  recorder  to  be  analyzed 
later.  Processing  of  the  data  was  done  on  a  dual  channel  400  B  omniferous  FFT  analyzer 
which  could  process  two  channels  of  data  simultaneously  in  real  time.  Excellent  resolution 
was  obtained  through  800-element  alias-free  frequency  functions  computed  from  dual  2048 
point  sampled  time  histories.  For  the  present  study,  100  ensemble  averages  of  each  function 
over  a  1  KHz  bandwidth  was  found  to  adequately  represent  the  time  independence  of  the 
stochastic  processes  being  investigated. 

The  results  for  the  pipe  and  annular  configurations  are  presented  in  Figs.  3  and  4 
respectively  for  a  typical  sensor.  The  non-dimensionalization  used  is  consistent  with  other 
investigators,  with  u;  denoting  the  radian  frequency,  Sm  one  half  the  diameter  ((r0  -  r,)/ 2 
for  annular  flow)  and  Uav9  the  average  bulk  velocity  obtained  through  conservation  of 
mass.  PSD  refers  to  the  power  spectral  density  (Pa2 /Hz)  and  p  is  the  fluid  density.  As 
illustrated,  the  level  of  facility  related  sound  (obtained  by  cross-correlating  two  coplan ar 
sensors)  is  approximately  5-10  dB  below  the  level  of  turbulent  fluctuations  over  1  KHz. 
This  level  increases  somewhat  in  the  annular  configuration  with  slight  fluctuations  in  the 
turbulent  spectrum  T,  as  might  be  expected  when  changing  to  the  more  restricted  flow 
regime.  Figure  5  shows  that,  for  both  pipe  and  annular  flows,  the  pressure  spectrum  scales 
with  outer  variables  for  all  velocities  investigated. 

Finally,  a  typical  run  of  pipe  flow  data  is  compared  in  Fig.  6  to  similar  results  in  the 
literature  with  excellent  agreement.  It  should  be  noted  that  Bakewell’s  (2)  data  has  been 
corrected  for  high  frequency  cancellation  according  to  the  method  proposed  by  Corcos 
(17),  whereas  the  present  data  and  that  of  Clinch  (5)  have  not.  The  excellent  comparison 
of  our  data  to  Bakewell’s  suggests  that  the  transducers  used  in  this  study  combined  with 
the  technique  presented  here  have  produced  excellent  resolution  and  accuracy  and  hence 
the  Corcos  correction  was  felt  unnecessary. 


CONCLUSIONS 


High  quality  wall  pressure  fluctuation  experiments  have  been  conducted  within  two 
internal  flow  configurations  with  the  following  basic  results. 
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(1)  A  new  signal  processing  technique  has  been  presented  which  utilizes  the  various 
subtracted  signals  from  three  coplan ar  sensors.  The  method  was  shown  to  represent  a 
more  straightforward  approach  than  those  previously  reported  . 

(2)  By  cross-correlating  two  coplan  ar  transducers  which  are  spatially  separated  such 
that  zero  turbulent  correlation  is  achieved,  a  measurement  of  the  background  facility  gen¬ 
erated  noise  levels  may  be  recorded. 

(3)  The  method  has  been  shown  to  be  easily  applied  to  estimating  finite  correlations 
which  may  or  maynot  exist  between  two  closely  spaced  transducers. 

The  effectiveness  of  this  method  is  limited  to  acoustic  noise  which  propagates  in  a 
plane  wave.  In  order  to  determine  the  relative  merits  of  applying  this  technique  to  higher 
order  propagation,  an  error  analysis  would  have  to  be  conducted.  Such  an  endeavor  is 
planned  in  the  near  future. 
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Figure  2.  Blow  down  water  tunnel  in  the  Laboratory  of  Computational  Physi 
and  Fluid  Dynamics  at  the  Naval  Research  Laboratory. 
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Figure  4.  Annular  pressure  spectra  at  Uavi  equal  to  6.5  m/aec. 
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Figure  5.  Pressure  spectra  comparison  of  the  turbulent  boundary  layer 
fluctuations  at  typical  flow  velocities  in  both  geometries 


Turbulent  pressure  spectra  comparison  at  a  How  speed  ot 
approximately  11. 0  m/sec  for  pipe  flow,  (see  refs.  2  and  5 
for  Bakcwell  and  Clinch  respectively.) 


